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Abstract 

We investigated the anomalous peak resistivity below the onset T c in un- 
derdoped YBCO, reminiscent of that observed in ID wires of conventional 
superconductors. We performed measurements of the angular dependence of 
resistivity p{9) in a magnetic field and the temperature dependence of resistiv- 
ity p(T), which exhibit a peak for B || ab-planes. This peak in p{T) disappears 
for B || c-axis. The width of the corresponding maximum in p{9) at 9 = 0° (B || 
ab-planes) decreases with increasing c-axis component of the field {B sin 9). The 
maximum in p{9) and p{T) decreases with an increasing applied transport cur- 
rent. We analyzed the data using three different models of resistivity based on 
2D resistor array, flux motion, and thermally activated phase-slips. Numerical 
calculations suggest that in a filamentary underdoped system, the phase-slip 
events could produce the anomalous resistivity close to T c . 

(PACS numbers: 74.76.-w, 74.40.+k,74.72.Bk) 
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1 INTRODUCTION 



Observation of a large resistive peak in the temperature dependence of resistivity 
p(T) just below the onset T c has been reported in crystals of high T c superconductors 
(HTSC) like {Nd,Pr) x ^Ce Q . lh Cu0^ y [|, YBCO(123) @, and BSCCO(2212) 0. In 
all these cases, the magnitude of the resistive peak is higher than the resistivity at the 
onset T c . The peak shows an anomalous behavior in a magnetic field applied along 
the c-axis. Its magnitude decreases with an increasing field, and in high enough fields 
the peak is completely suppressed. The applied transport current has the same effect 
on the peak i.e the peak's magnitude decreases with an increasing applied current. 

These phenomena are of considerable interest because of their striking qualita- 
tive similarity to those observed in conventional superconductors (LTSC) like super- 
conducting mesoscopic Al wires (§ §, thin films of Al g 0|, (NbV)N, NbV, VN, 
(NbTi)N ||, and disordered metallic glasses of Zr 60 Cu 60 0. Mosqueira et al re- 
ported the resistive anomalous peak in p a b(T) of YBCO crystals of T C (R = 0) = 89K. 
This anomaly was eliminated by successive annealing of the sample in oxygen. This 
annealing also led to an increase of T C (R = 0) up to 90.3 K (close to the opti- 
mal doping). The authors concluded that the peak could be related to very small 
T c -inhomogeneities non-uniformly distributed in the crystal. They performed com- 
puter simulations of the temperature dependence of the anomalous resistivity using 
the model of two-dimensional electric circuit: an array of resistors whose resistivity 
depends on temperature. The non-uniformly distributed T c -inhomogeneities were in- 
troduced by assuming that these resistors have different (higher or lower) T c . The 
authors stated that the uniformly distributed T c -inhomogeneities at large length scales 
broaden the resistive transition only, and do not produce a peak. Similar approach 
was introduced earlier by Vaglio et al |§ to explain the resistance-peak anomaly 
in non- homogeneous thin films of (NbV)N, NbN, VN, and (NbTi)N. They concluded 
that "the current redistribution" caused by the sample's inhomogeneity, is responsible 
for the observed phenomena. 

Current redistribution effects in an inhomogeneous sample were also considered 
by Nordstrom and Rapp |J in their interpretation of the resistive peak anomaly 
in superconducting amorphous thick films of Zr eQ Cu i0 . Kwong et al observed 
an anomalous peak in the resistive transition of a 2D 25 nm thick aluminium film 
containing regions of different but comparable, transition temperatures. Disordered 
regions of lower T c were produced by the reactive-ion etching process. Their observa- 
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tion seems to support earlier interpretations based on T c -inhomogeneities and current 
redistribution effects. The authors stated, however, that the anomaly could originate 
from a discontinuity of the superconducting potential at the normal-superconducting 
metal (N-S) interface, and for superconducting electrodes placed sufficiently close to 
the interface, this potential exceeds the normal-state value. Spahn and Keck found 
that the anomaly appears in 2D Al films with thickness between 13 and 40 nm. They 
argued that this effect could be caused by an interaction between the superconducting 
fluctuations and the conduction electrons. 

Extensive studies of the resistive anomaly were also performed on ID Al strips 
with a width less than the coherence length and the magnetic penetration depth, by 
Santhanam et al || and Moshchalkov et al |J. Santhanam et al argued that the Al 
wire could be treated (at temperatures close to T c ) as a coherent region comprising 
normal (N) and superconducting (S) phases. The resulting N-S interface gives rise to 
a quasiparticle charge imbalance induced by the bias current, and consequently to the 
observed changes in resistivity. Moshchalkov et al performed quantitative analysis of 
the anomaly using Langer-Ambegaokar (LA) [|l(J and McCumber-Halperin (MH) | 1 
models of the thermally activated phase-slips of the superconducting order parameter. 
LA-MH models were adopted with the modification which assumes that in quasi-lD 
superconducting wires the normal current and the supercurrent can only flow in series, 
and the total resistance is the sum of the normal resistance Rn and the phase-slip 
resistance Rs- Good quantitative agreement between the experimental data and the 
calculated resistance R(T) was obtained. 

Crusellas et al and Han et al || proposed that the anomaly in p a b(T) of 
(Pr, Nd)i .^CeQ ^CuOi-y and BSCCO(2212) crystals is the manifestation of a quasireen- 
trant behavior, which results from the intrinsic granularity. Han et al rejected 
the explanation based on non- uniformly distributed T c -inhomogeneities (Ref. 0), be- 
cause of the observation of an anomalous peak in the I-V characteristics, which 
were measured at different magnetic fields. However, Crusellas et al |I[ stated that 
the anomaly is strongly influenced by the distribution of defects, after it was ob- 
served that high temperature annealing reduces the size of the resistive anomaly in 
(Pr, Nd)i 85 Ceo,i 5 Cu0^ y crystals. 

Briefly, the interpretation of the anomalous resistive peaks in HTSC concentrates 
on two possible sources of this effect: non- uniformly distributed T c -inhomogeneities 
and intrinsic granularity. The explanation of the anomalous resistivity in LTSC films 
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took into account the effects of T c -inhomogeneities (and related current redistribu- 
tion), N-S interfaces, and the interaction between superconducting fluctuations and 
the conduction electrons. It was also suggested that the anomaly in ID LTSC (Al) 
strips (wires) originates from the presence of N-S interfaces and/or thermally acti- 
vated phase-slips of the order parameter. 

These various interpretations are the source of a number of unanswered questions: 



1. According to Browning et al [0 in YBCO single crystals of T c =93 K and tran- 
sition width of AT C =0.2 K, large variation in the oxygen content 7 — 5 can 
occur across the sample as revealed by high resolution scanning x-ray diffrac- 
tometry (which was performed using 10/im wide x-ray beam). 7 — 5 in these 
crystals ranges between 6.80 and 7.00, which corresponds to a change of T c by 
about 10 K. In spite of these non-uniform T c -inhomogeneities, the crystals have 
small resistivities (p ~ 40pQcm at 100 K) and do not show any resistive peak 
anomalies at the onset T c in p(T). These results throw doubt on whether the 
2D resistive model alone (as proposed in Ref.0 for YBCO) can explain the 
observed anomalies. 

2. The resistive anomalies observed in LTSC films and wires (strips) are similar and 
their interpretation suggest the link between the presence of inhomogeneities 
(T c -inhomogeneities, N-S interfaces) and the superconducting fluctuations, in- 
cluding phase-slips of the order parameter. Could this explanation be also 
applied to HTSC? 

3. Moshchalkov et al || introduced the phase slip resistivity ( according to the ID 
LA-MH model |T0| , [TT|| ) combined with the normal state resistivity in order to 
explain the anomalous resistive peak in ID aluminum wires. Experiments by 
Browning et al [see (1)] suggest filamentary phase separation and filamentary 
flow of the current in some YBCO crystals with sharp superconducting transi- 
tions, which do not show resistive anomalies. Does this mean using the anal- 
ogy to LTSC that the presence of non-uniformly distributed inhomogeneities 
in HTSC is the necessary but not sufficient condition to observe the resistive 
anomaly? What is the other condition? Could this be a ID current flow in an 
inhomogeneous system? 

4. What is the contribution of the magnetic flux motion (pinning) to the observed 
resistive anomalies in HTSC? 
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In order to answer these questions new experiments are needed. We decided to 
perform measurements of the angular dependence of resistivity in a magnetic field. 
This decision was stimulated by the experiments done on Pri^Ceo^Cu^-y crystals 
IJ, in which the resistive anomaly was investigated for two different directions of 
the applied magnetic field, namely along the c-axis and along the ab-planes. The 
effect of the magnetic field on the anomaly was completely different for these two 
orientations. Taking into account the fact that the presence of the non-uniform T c - 
inhomogeneities is the necessary but not sufficient condition to observe the resistive 
anomalies, we investigated the temperature dependence of resistivity on a large num- 
ber of YBCO samples. We decided to study YBCO thin films, both optimally doped 
and underdoped, because they are readily accessible from different research groups 
and can be deposited on various substrates using several different deposition tech- 
niques. The bridges for the resistive measurements can be made relatively easily on 
thin films using standard photolithographic techniques. The resistive peak anomaly 
and the reduction of its magnitude with an increasing magnetic field and an increas- 
ing applied current, were observed in an underdoped film after investigation of fifteen 
YBCO films of T c ranging between 79 and 90.5 K. This film was then used to perform 
detailed measurements of the angular dependence of resistivity in a magnetic field. 
The resulting experimental data were analyzed using three different models: two di- 
mensional resistor model, magnetic flux motion model, and the LA-MH thermally 
activated phase -slip theory. 

2 EXPERIMENTAL PROCEDURE 

2.1 Sample Preparation 

C-axis oriented YBCO thin films were prepared using off-axis rf magnetron sputtering 
and laser ablation from stoichiometric YBa2Cus07-s targets of 99.999% purity. Films 
were deposited on three different types of substrates: SrTiO^, LaAlO^, and sapphire 
(with Ce0 2 buffer layer). 

We investigated fifteen YBCO thin films (both underdoped and close to the opti- 
mal doping) of various zero-resistance transition temperatures (between 79 and 90. 5K) 
and thicknesses (between 100 and 600nm). YBCO films were patterned, using con- 
ventional photolithography and wet etching technique, into a form of a 30 — 60fim 
wide and 6.4 mm long strips with six measurement probes. Large area silver contacts 
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were deposited on the film by rf magnetron sputtering in order to minimize Joule 
heating. Copper leads were attached to silver contacts using mechanically pressed 
indium. The distance between voltage probes was 0.4 mm. 

The anomalous resistivity was observed in an underdoped YBCO thin film. This 
film (140 nm thick), was deposited on a (1000) oriented sapphire substrate (with CeOi 
buffer layer) using laser ablation technique. The sample exhibits a vanishing zero-field 
resistivity at T c = 81.7K and has a room temperature resistivity p^qqk = 34.2pQcm 



' P300K 



2.4). 

PlOOif ' 

X-ray diffraction (XRD) data of this film showed the pattern of a characteristic 
stoichiometric c-axis oriented YBCO film. The data did not reveal any impurity 
phases. The XRD data gave a c-axis lattice spacing of 11.70 A, which corresponds to 



an oxygen content of about 6.8 and T C (R = 0) of about 80 K 13 



2.2 Measurement Procedure 

The investigation of the resistive anomalies was based on the following measurements: 
(a) the measurement of the temperature dependence of resistivity p(T) between room 
temperature and T C (R = 0) in zero magnetic field; (b) the measurement of the tem- 
perature dependence of resistivity p(T) between the onset T c and T C (R = 0) in an 
external magnetic field applied either parallel or perpendicular to the ab-planes; (c) 
the measurement of the angular dependence of resistivity p(6) as a function of the 
angle 9 between the ab-planes and the direction of the fixed applied magnetic field 
at fixed temperatures between the onset T c and T C (R = 0); (d) the measurement of 
p{6) as a function of the magnitude of the magnetic field B and the applied transport 
current density J. The angular measurements were performed by rotating a copper 
sample holder about its vertical axis in a horizontal magnetic field up to 1 Tesla, 
using a combination of a step-motor and backlash-free gear reducer. The angle was 
accurately monitored by an 8000-line optical encoder attached to the sample, whose 
angular resolution was 0.045°. The film was mounted with the c-axis perpendicu- 
lar to the sample holder's vertical axis of rotation, which allowed one to change the 
magnetic field direction in a plane parallel to the c-axis. 

Resistivity was measured using the standard dc four-probe method. The current 
was applied to the sample in the form of short pulses (of duration less than 200 
ms) in order to reduce Joule heating. A dc current reversal was used to eliminate 
the thermal emf in the leads. The voltage was measured using a Keithley 2182 
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Nanovoltmeter in tandem with a Keithley 236 Current Source, with the nanovoltmeter 
working as the triggering unit. The nanovoltmeter was operated in a mode (known as 
Delta mode) which allows the measurement and calculation of the voltage from two 
voltage measurements for two opposite directions of the current. Temperature was 
monitored by a carbon-glass resistance thermometer and an inductanceless heater and 
was controlled to better than HOmK for each single angular sweep in a magnetic field. 
This was achieved by rotating the sample very slowly in the magnetic field in order 
to reduce variations in the emf in the heater which could disturb the temperature 
reading. The term "resistivity" is used in this paper to denote the quantity E/J 
(where E is the electric field and J is the transport current density), and it does not 
imply an ohmic response. 

All measurements were carried out with the transport current J parallel to the ab- 
planes for two different orientations of the magnetic field with respect to the current. 
For the first one, the field was rotated in a plane perpendicular to the current direction 
while for the other one the field was rotated in a plane parallel to the current and the 
c-axis directions [see Fig.|l|(a)]. All measurements were done in a field cooling (FC) 
regime, with the magnetic field applied to the sample at a temperature above the onset 
T c , followed by a slow cooling down to the required temperature of measurement. 

3 EXPERIMENTAL RESULTS 

3.1 Temperature dependence of resistivity 

The temperature dependence of resistivity p was measured over a temperature range 
of 78-300 K in a zero magnetic field. For a temperature range (78-90 K) close to 
T c , p was recorded for different orientations of the magnetic field with respect to 
the direction of the current density J. Fig.[I|(a) shows two possible orientations of 
the magnetic field B with respect to J and the ab-plane of the film. The figure on 
the left illustrates the case in which the field B was rotated in a plane parallel to 
the direction of J while the one on the right represents the case in which B was 
rotated in a plane perpendicular to J. For both orientations B was rotated in a 
plane parallel to the c-axis. The first configuration is denoted as B || J and the 
second one as B _L J. Fig.|l](b) shows the temperature dependence of resistivity p(T) 
for a temperature range of 79-88 K measured in a zero field and in 0.68 T. The 
measurements of p(T) in the field were carried out for the following orientations of 
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B with respect to J: B || J and B _L J with B parallel to the ab-plane (9 = 0°), and 
for B _L J with B parallel to the c-axis (9 = 90°). The onset transition temperature 
(onset T c ) is defined as the temperature above which the resistivity does not respond 
to the change in both magnitude and direction of the magnetic field [see Fig.|l](b)]. 
p(T) below the onset T c could be divided into three regions. Each region is identified 
according to the response of p(T) to the change in the direction of the magnetic 
field from the ab-planes (9 = 0°) to the c-axis (9 = 90°). Region II represents a 
temperature range between 82.8 K and 83.5 K over which p(T) exhibits a peak (of 
magnitude larger than that of p(T) at the onset T c ) in a zero magnetic field, and for 
B || ab-planes with B || J and B _L J orientations. Note a clear separation between 
the peak and the onset T c . In this region, behavior of p(T) changes dramatically 
upon rotating the field from the ab-planes to the c-axis. In regions I and III, p(T) 
was observed to increase when B is parallel to the c-axis, while in region II (the peak 
region) p(T) is completely suppressed by the magnetic field B || c-axis. For B | 
ab-planes, p(T) in regions I and III is independent of the magnitude of B and the 
angle between B and J. However, in region II, p(T) is independent of the magnitude 
of B only for B || J orientation. In this region, p(T) was found to decrease with an 
increasing applied current density J. The temperature dependence of resistivity for 
B || c-axis was measured in different magnetic fields. For fields above 0.1 T, the peak 
in region II is completely suppressed. The resistivity between the onset T c and the 
room temperature exhibits a linear temperature dependence. 

3.2 Angular dependence of resistivity: Effect of temperature 
and magnetic field 

The angular dependence of resistivity p{9) was measured in a constant magnetic 
field at different temperatures in regions I, II, and III, for both B || J and B _L J 
orientations for the angular range from —20° to +20°. The measurements revealed 
minima in p(9) at 9 = 0° in region I and III, and a maximum in region II (see Fig.|2|). 
Fig.||](a) shows p{9) in region I as a function of temperature between 82.42 K and 
82.82 K in a magnetic field of 0.68 T. At a temperature of approximately 82.82 K, 
which corresponds to the border line between region I and II in Fig.[l|(b), there is a 
crossover from a minimum to a peak in p{9). This peak grows with an increasing 
temperature reaching a maximum value at 83.23 K (see Fig.|2|(b)). 

The second crossover from a maximum to minimum can be seen at 83.50 K, 
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Figure 1: (a) Two configurations of B with respect J that were used during the 
measurements of the angular dependence of resistivity p{6) in a magnetic field: B 
is rotated in a plane parallel to both J and the c-axis (left side), or B is rotated in 
a plane parallel to the c-axis but perpendicular to J (right side), (b) Temperature 
dependence of resistivity for YBCO thin film measured in a zero and 0.68 T fields at 
different orientations. Regions I and III denote the temperature ranges over which 
p(T) is independent of the magnitude of B and the angle between B and J for B || 
ab-planes. In these regions p(8) displays a minimum at 9 = 0° (B parallel to the 
ab-planes) [see Fig.^|(a) and (c)]. Region II represents the temperature range over 
which p(T) exhibits a peak (of magnitude larger than p(T) at the onset T c ), for B=0 
and for both B || J (6 = 0°) and B _L J (0 = 0°) orientations. Rotating the field from 
the ab-planes (6 = 0°) towards the c-axis (6 = 90°), leads to an increase in p(T) in 
regions I and III and to a suppression of the peak in region II. 
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Figure 2: p(6) measured in 0.68 T for a temperature range between 82.52 K and 83.83 
K spanning the three regions I, II, and III. Note the change in p(9) at 9 = 0° from a 
minimum in region I (a) to a maximum in region II (b) and then back to a minimum 
in region III (c). In regions I and III, identical behavior of p{9) has been observed for 
B || J and B _L J orientations, whereas in region II the magnitude of p(9) depends 
on those orientations. 
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which corresponds to the border line between regions II and III in Fig.[T](b) . We have 
measured p{9) for both B _L J and B || J in all three regions. While in regions I 
and III the minimum in p{9) is independent of the orientation of B with respect to J 
(i.e for B _L J and B || J), in region II the magnitude of the peak depends on these 
orientations and p(9) B ±j > p(9)b\\j- 

The measurements of p(6) over an angular range between —30° and 210° revealed 
sharp maxima for B || ab-planes (9 = 0° and 9 = 180°) and a smaller broad maximum 
for B || c-axis (9 = 90°) (see Fig.|). p(9) at 9 = 0° and 9 = 180° is about 30% larger 
than that for 9 = 90°. Moreover p(9) has minima at 9 = 35° an 9 = 145° for all fields. 




B // ab-plane B // ab-plane - 



-30 30 60 90 120 150 180 210 
e (degrees) 



Figure 3: p(9) measured in 0.68 T at 83.23 K for an angular range —30° < 9 < 210°. 
Note the peaks at 9 = 0° and 9 = 180° (B || ab-planes) which are approximately 30% 
higher than the maximum at 9 = 90° (B || c-axis). 

The angular dependence of resistivity p{9) was also measured at a constant tem- 
perature in different magnetic field in regions I, II, and III, for both B || J and B _L J 
orientations and for the angular range between —20° and +20°. Fig.f| shows the an- 
gular dependence of resistivity, at a temperature of 81.43 K (region I), measured in 
different magnetic fields for both B || J and B _L J orientations. The data for these 
two orientations are identical which implies that p(9) is independent of the angle 
between B and J. p(9) at 9 = 0° is almost independent of the magnitude of B. The 
width of this minimum [defined as half width at half minimum (HWHM)] decreases 
from HWHM=2.3° in a field of 0.17 T down to HWHM=2.0° in 0.86 T. The depth 
of the minimum increases with an increasing field. The results of the measurements 
of p(9) in region III is identical in all aspects to those obtained in region I. 
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Figure 4: (a) Angular dependence of resistivity measured at a temperature of 81.43 
K (region I) in different fields for B close to the ab-planes. The data reveal a minimum 
in p{6) at 9 = 0° (B parallel to the ab-planes). The width of the minimum decreases 
with an in increasing field. An identical behavior has been observed for both B _L J 
and B || J orientations. 

Fig.[| presents p(9) measured at a temperature of 83.13 K (in region II) in different 
magnetic fields for both B _L J and B || J orientations. The width of the peak in 
p{6) decreases with an increasing B for both B _L J and B || J orientations. The 
magnitude of the peak in p{9) for B || J is almost independent of the magnitude of 
B, however it increases with B for B _L J. A decrease of the peak's width with an 
increasing magnetic field means that within a certain angular range (| 9 |> 1.5° for 
B _L J and | 9 |> 0.3° for B || J), p{9) decreases with an increasing field. 

3.3 Angular dependence of resistivity: Effect of the applied 
current 

The angular dependence of resistivity was measured also as a function of the applied 
current density at a constant temperature and magnetic field. Fig.^| presents the 
measurements of p(9) in region II for a wide range of applied current density J, 
from 0.9 to 69.4 kA/cm 2 , in a field of 0.68 T and at a temperature of 83.03 K. For 
angles \ 9 \> 2°, p(9) increases non-linearly with an increasing J, but for small angles 
| 9 \< 2° it decreases with an increasing J. In the angular region for | 9 |< 1°, starting 
at small current density (J < 11.6kA/cm 2 ), the peak height initially decreases with an 
increasing current, but for J larger than 23.1k Aj cm 2 , a minimum in p(9) develops. 
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Figure 5: Angular dependence of resistivity p(9) for —5° < 9 < 5°, measured in 
different applied magnetic fields at a fixed temperature of 83.13 K (region II). (a) 
p(9) for B _L J orientation, where p(9) increases with an increasing B for 9 < 1° 
but decreases with an increasing B for 9 > 1°. (b) p(9) for B || J orientation; p at 
9 = 0° is almost independent of B. Note that p(9) decreases with an increasing B for 
9 > 0.5°. 

The dependence of p{9) on J is essentially the same for both B || J and B _L J 
orientations. 

Fig.fj] shows p{9) measured in region I for a wide range of J in a field of 0.68 T 
and at a temperature of 81.83 K. Effect of the current on the minimum is different 
from that observed in region II. The minimum at 9 = 0° decreases with an increasing 
J. The dependence of p{9) on J is identical for both B || J and B _L J orientations. 
Similar dependence of p(9) on J was observed over the temperature range in region 
III. 



4 DISCUSSION 

The experimental data for p(T, B) obtained for the underdoped YBCO film are qual- 
itatively similar to those observed before in HTSC [1], 0, The anomalous resistive 
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Figure 6: Angular dependence of resistivity p(9) measured as a function of applied 
current density J for a field of 0.68 T and a temperature of 83.03 K (region II) . For 
angles | 9 \> 2° , p(9) increases with an increasing J, but for small angles | 9 |< 2° 
the opposite happens, where a minimum starts to develop with its width increasing 
with an increasing applied current. 

peak is located at a temperature approximately 1.8K lower than the onset T c . The 
peak disappears when a magnetic field is applied along the c-axis of the film. Also the 
magnitude of the peak decreases and its position shifts to lower temperatures with 
an increasing applied transport current. Our measurements of the angular depen- 
dence of resistivity in a magnetic field provided very valuable additional information, 
which allows us to understand better the physics of the anomalous resistivity. The 
measurements of p(9) in a magnetic field as a function of temperature revealed sharp 
minima in resistivity at 9 = 0° (B parallel to the ab-planes) at temperatures below 
and above the resistive peak in p(T), and a sharp maximum at 9 = 0° at the peak's 
temperature (see Fig.|2|). 

The data for p(T) and p(9) were used to distinguish between different interpreta- 
tions of the resistive anomaly. We considered two-dimensional resistor model, mag- 
netic flux motion, and thermally activated phase-slips. 
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Figure 7: (a) Angular dependence of resistivity p(8) measured as a function of an 
applied current density J in a field of 0.68 T and at a temperature of 81.83 K (region 
I). The sharpness of the minimum at 9 = 0° decreases gradually with an increasing 
J. (b) Expanded view of p(8) measured at J = 57.9kA/cm 2 . 

4.1 Two-dimensional resistor model 

In an inhomogeneous superconductor different parts of the film or the crystal can 
have slightly different transition temperatures. In this case the superconductor may 
be modelled as an electrical circuit- array of different resistors. The anomalous peak 
in p(T) is produced by solving numerically, through the standard matrix method, the 
electrical circuit equations. This model was used to analyze the data, taken in a zero 
magnetic field, for anomalous resistivity in p a b{T) of YBCO crystal by Mosqueira et 
al [Q], and in p(T) of Nb-based LTSC films by Vaglio et al ||. The calculated p(T) 
agrees with the experimental data. 

Mosqueira et al also attempted to explain the suppression of the anomalous 
peak by a magnetic field in YBCO crystals using the resistor model. They argued that 
the reduction of the peak and its shift to low temperature is caused by the broadening 
of the superconducting transition AT C (B*) in a magnetic field B* (and not just by the 
shift in T c ). The peak disappears when the superconducting transition broadens by 
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Table 1: Ap/p B , AT C (B*), and the magnetic field B* at which the peak disappears, 
for different superconductors. 



Material 


Ap/pB 


B*(T) 


AT C (B*)(K) 


Ref 


YBCO crystal 


0.29 


0.3 


2 


§ 


YBCO film 


0.44 


0.08 


0.5 


this work 


BSCCO crystal 


1.56 


0.01 


~ 2 


§ 


NdCeCuO crystal 


0.28 


0.1 


r*j 4 


E 


PrCeCuO crystal 


0.16 


0.7 


< 0.5 


E 


Al wires 


0.17-0.56 


0.001 





E 



about 2K in a field of 0.3T. According to the resistive model, the broadening required 
for the peak to be eliminated from p(T) should be material-independent provided 
that the ratio Ap/p B = (p P — p B )/ p B (where pp is the peak's resistivity and p B is the 
resistivity measured in a magnetic field B at the peak's temperature in the absence 
of the peak) does not change. The experimental data in Refs |jj |2|, [3], |[] show that 
this is not the case. Table [3] lists the data for Ap/p B , AT C (B*), and the magnetic 
field B* at which the peak disappears, for different superconducting materials. 

These data reveal that the disappearance of the peak from p(T) in a magnetic 
field is not related to the broadening AT C (B*) of the superconducting transition. Our 
data obtained on YBCO film (see Table [!]) support this conclusion. 

4.2 Flux-motion 

In order to find the contribution of the magnetic flux motion to the observed resistive 
anomaly, we performed the measurements of p{T) and p{9) for two different orienta- 
tions of the current relative to the magnetic field i.e for B _L J and B || J orientations 
(see Fig.|l|). The magnitude of the peak in p(T) (region II in Fig.[l]) increases when 
the field is parallel to the ab-planes and perpendicular to the current i.e for B _L J 
(compared to the case for B || J). This situation corresponds to the maximum Lorenz 
force acting on the flux lines along the ab-planes. The angular dependence of p in a 
magnetic field (Fig.0) reveals a maximum in region II, but sharp minima at temper- 
atures below and above the peak (regions I and III in Fig.[l]). A very sharp minimum 
in p{9) at 9 = 0° (B || ab-planes) was seen previously in a YBCO single crystal 
by Kwok et al and interpreted as due to the lock-in transition of the flux lines 



16 



trapped between the planes. For a system of weakly coupled C11O2 layers one expects 
a maximum resistive dissipation for B || c-axis and a minimum for B || ab-planes. An 
increase in p{9) when the field is rotated from the ab-planes to the c-axis is normally 
attributed to the intrinsic anisotropy of the material. We found that the minima in 
p(9) at 9 = 0° (regions I and III) are independent of the orientation of the current 
relative to B (see Fig.^), suggesting very strong flux lock- in mechanism when B is 
parallel to the ab-planes. For 9 > 0° resistivity could arise via the nucleation and 



motion of kinks along the vortex lines ||14|| . In this case, one could also describe the 
tilted vortex line as the combination of Josephson strings aligned along the ab-planes 
and mobile pancakes (vortex segments along the c-axis). If the coupling between the 
pancake vortices is weak, the Lorenz force acting on these vortices and consequently 
their motion, should be independent of the direction of the transport current in the 
ab-planes. The measurement of the minimum in p(9) also revealed an increase of 
the resistivity with an increasing transport current in the ab-planes (Fig.[7]), which is 
independent of the orientation of the current relative to B. This result suggests that 
the motion of the pancake-vortices in the ab-planes is responsible for the observed 
increase of p(6) for 9 > 0° in regions I and III. The maximum in p(6) at 9 = 0° at 
temperatures corresponding to region II in p(T) (Fig.|5|) depends on the orientation 
of the current relative to B. For B _L J orientation, the maximum is higher than 
that measured for B || J orientation. This behavior is different from that observed in 
regions I and III and therefore it provides additional argument that the peak in p(T) 
can not be explained by the 2D-resistor model alone. It also suggests that the un- 
known dissipation in region II weakens flux lock-in between the planes. Subtracting 
the maximum in p(6) at 9 = 0° for B || J from that measured for B _L J (see Fig.[8p 
gives p(9) with a minimum similar to those observed in regions I and III, which are 
caused by flux motion. 

The measurement of the maxima in p(9) at 9 = 0° as a function of magnetic 
field for B || J and B _L J orientations shows that the maxima become sharper 
(i.e their width decreases) with an increasing magnetic field. For both B || J and 
B _L J orientations, and for 9 > 1.5°, the resistivity at a fixed 9 decreases with an 
increasing field (Fig||). On the other hand, the maximum in p{9) at 9 = 0° also 
decreases with an increasing transport current (Fig.|6]). This reduction in resistivity 
can not be explained by the flux motion. Chaparala et al |I3fl observed a small 
maximum in p(9), when the magnetic field was oriented parallel to the ab-planes in 
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Tl (2212, 1223) and Bi(2212) crystals. The authors did not present any data for the 
corresponding temperature dependence of resistivity. The maximum in resistivity 
was attributed to the formation and motion of the c-axis-oriented vortex-antivortex 
segments of the flux lines parallel to the ab-planes. They assumed that the maximum 
is created as a result of the interplay between the density n s and the velocity v s of 
the vortex-antivortex segments. The resistive potential difference V is proportional 
to the product of these quantities. According to the experimental observation V oc 



(n s v s )g = Qo is larger than V oc (n s v s )o >0 o. Chaparala et al |15| argued that at 9 = 0°, 
in spite of the small density of the vortices, n s v s is large because of the high velocity 
of the newly created vortex-antivortex pairs. At 9 > 0°, n s is large but v s is small, 
so (n s v s )g =0 o > (n s v s ) e>0 o. According to this interpretation, increasing the applied 
transport current should increase the Lorenz force on these pairs and consequently 
increase their velocity. This leads to an increase in the resistive dissipation and to 
the growth of the maximum in p(9) at 9 = 0° with an increasing current. Our data 
revealed a reduction of the maximum in p{9) at 9 = 0° with an increasing current (see 
Fig.p), which eliminates the vortex-antivortex model as a possible explanation of the 
resistive anomaly. The absolute values of the resistivity in the peak observed on p(T) 
curve is higher than the resistivity at the onset T c (85 K), defined as the temperature 
above which the resistivity is independent of the magnitude and direction of the 
applied magnetic field (see Fig.|I|). The resistive dissipation due to a vortex motion 
can reduce the critical current density to zero, reaching the normal state resistivity, 
but it can not exceed this value. 



4.3 Phase-slip model 

Discussion of the resistive-peak anomaly in the previous sections indicates that 2D- 
resistor and flux motion models alone can not fully account for the origin of this phe- 
nomenon. Regarding LTSC, Moshchalkov et al || argued that the resistive anomaly, 
seen in ID Al wires, originate from thermally activated phase slips, and the observed 
resistive peak at the onset T c is the result of the phase-slip resistivity and the normal 
state resistivity acting in series. Observation of the similar resistive peak anomaly in 
LTSC disordered films implies that in some disordered systems filamentary flow of the 
transport current could occur through ID constrictions (channels). We believe that 
this could also happen in HTSC samples. Browning et al |12|] revealed that in spite 
of a large variation of the oxygen content (7 — 5 = 6.8 — 7.0) measured across YBCO 
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Figure 8: (a) Comparison between p{9) measured for B _L J and B || J orientations 
in 0.68T at 83.13K. (b) The difference Ap(6) between the peaks in p{9) for B _L J 
and B II J orientations measured in 0.34 T and 0.68 T. 



crystals, they still display sharp superconducting transitions (~ 0.2/T), high T c , and 
low resistivity. This implies filamentary flow of the current in the samples. However, 
the resistive peak anomaly is absent in the samples which could mean that the fila- 
mentary flow alone is not sufficient to produce the resistive anomalies. We conclude 
by analogy to the case of LTSC disordered films that the thermally activated phase- 
slips could produce such anomaly if the filamentary flow of the current occurs through 
ID constrictions. This could happen more likely in underdoped HTSC samples due 
to phase separation-induced disorder. It should be also noted that our data show all 



qualitative basic characteristics expected by the LA-MH phase slip model ITT], ITU 



According to this model phase slip events lead to the appearance of a resistance in 
ID superconducting wires below T c . During a phase slip event, thermal fluctuations 
reduce the superconducting order parameter, defined as ip(x) =\ ip(x) \ e l ^ x \ where 
<f)(x) is the phase, to zero at some point along the wire momentarily disconnecting the 
phase coherence. This allows the relative phase across the wire to slip by 2-ir (before 
4>(x) recovers its finite value), resulting in a resistive voltage. 

For a ID thin wire with a transverse dimension d £ and d <C A, the LA-MH 
theory predicts that the appearance of a resistance in the superconducting state is 
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mainly determined by thermally activated phase slips events as the system is passing 
over a free energy barrier AF (the difference in free energy between the normal and 
superconducting states) proportional to the cross-sectional area A of the wire; 

AF = ^[^(T)tf c 2 (T)/8vr], (1) 

where H C (T) is the thermodynamic critical field. 

In the absence of the current, phase slips by ±2tt are equally likely, and this 
results in a fluctuating noise voltage with a zero net dc component. The result of 
the application of a current to the wire is to make the phase jumps more probable in 
one direction than in the other. The different jump rates arise from a difference SF 
in the energy barrier for jumps in two directions and this difference stems from the 
electric work / IV dt done in the process. For a phase slip of 2ir, the energy difference 
is 5F = AF + — AF_ = ^I s = 4>qI S) where (fi = ^- is the superconducting flux 
quantum. 5F = <p I s should be larger than the thermal energy k B T, which defines 
the characteristic current I x = above which most phase slips go in the driven 
direction and the resistance is nonlinear [[TIJ. I\ sets a lower limit on the applied 
current I s . The upper limit is set by the critical current I c which is the mean- field 
critical current given by I c = ^-^/f^S and 

I 1 = ^L<I S< I C (2) 

The average voltage V s arising from the phase slip events is determined by the 
number of these events in the sample [N(T) = L/£(T), where L is the length of 
the wire], a characteristic time r(T), Boltzmann factor exp(— AF(T)/kBT), and the 
factor sinh(I s (f) /2k B T) derived from the difference SF in the energy barrier for the 
+27r and — 2n phase jumps. V s is determined by 

AF(T)~ 



V s = 2<j> Q.{I B , T) exp 



k B T 



smh ( ^ ) (3) 



where AF(T) = AF (T) + (iy/ 2 I*k B T '/torlJc and Q{I S ,T) is an attempt frequency 
which can be approximated as 



n(/.,T) = ^J^(i-^ ( 4 ) 

V s ' ; r(T) V k B T y 3// V ; 

It is very important to emphasize the fact that the energy being supplied during 
the occurrence of these phase slips at a rate of IV, is dissipated as heat rather than 
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converted into kinetic energy of supercurrent, which would otherwise soon exceed the 
condensation energy JTo| . 

The magnetic field dependence of the phase-slip event does not appear explicitly 
in Eq.(|), however its effect on the phase-slip voltage appears through the dependence 
of the critical current I c on B [I c oc (1/B)]. Phase-slip resistivity is present over a 
range of the applied current I s between I x and I c according to Eq.(|2]). The applied 
current I s per ID current channel should be larger than Ii = ksT/(j) . If I s is too close 
to I c , the phase-slip events are less likely to occur. Also, reducing I c while keeping I s 
fixed leads to the reduction of the phase-slip events and consequently the voltage V s . 
For an anisotropic superconductor, increasing the magnitude of the c-axis component 
of B by increasing the angle 9 and/or the magnitude of B, reduces I c . 

The angular dependence of resistivity in a magnetic field measured over a temper- 
ature range between 82K and 84K (Fig||) points out different origin of resistivity in 
the peak in p(T) (Fig.|l|), in comparison to that at temperatures below and above the 
peak. Therefore p{T) could be treated as a superposition of the peak and the normal 
resistivity near the transition, which increases almost linearly with temperature be- 
tween T C (R = 0) ~ 82K and the onset T c ~ 85K. We considered the possibility that 
the resistive peak originates from thermally activated phase-slips and attempted to 
perform numerical calculations of the phase-slip resistivity using the modified LA-MH 
theory. We assumed that in an underdoped HTSC sample the current flows through 
n parallel superconducting filaments of length 0.4 mm (which is the distance between 
the voltage contacts). The width w (in the ab-planes) and the thickness t (along 
the c-axis) of a filament were chosen to be 2.0 nm and 1.0 nm, respectively. These 
values are much smaller than the coherence length in the ab-planes and along the 
c-axis at temperatures close to T c , and therefore the filaments can be treated as ID 
wires. In the system of n parallel superconducting filaments, one could expect that 
the phase-slip event occurring in a filament would affect the superconducting state 
of the neighboring filaments, because the coherence length at temperatures close to 
T c is much larger than the spacing between the filaments. On the other hand, in an 
underdoped system, one could also expect that along each filament the superconduct- 
ing regions are interrupted by segments of normal resistance Rn, so that the total 
resistance of the filament Rf is the sum of the phase-slip resistance and the normal 
state resistance acting in series: Rf = n s R s + Rn, where n s is the number of the su- 
perconducting segments of an average length l s and an average phase-slip resistance 
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R s . The arguments presented above suggest that for a system of n parallel filaments, 
the formula for the condensation energy AF (Eq. ([!])) for phase-slip events in a ID 
wire should be modified to reflect the phase slip events in the whole system of n ■ n s 
segments. We assumed AF in the form: 

AF = ^[AU(T)H 2 c (T)/8n] ■ n ■ n s , (5) 

where n is the total number of filaments, and n s is the average number of supercon- 
ducting segments per filaments (the average number of phase-slip centers per fila- 
ment). The Ginzburg-Landau expression for H C {T) = H C (0)(1 — T/T c ) = if c2 (0)(l — 
T/T c ) I \/2k, and £, a b(T) = £ a fc(0)/(l — T/Tc) 1 / 2 close to T c were used. The phase- 
slip voltage V s was calculated using Eq.@ and the modified attempt frequency 
Q(I S ,T). The number of the phase-slip events in a segment was given by N(T) = 
ls/£,ab{T)r(T) = •jL/^ a b{T)T{T), where 7 = l s /L. 7 and n s were treated as the fitting 
parameters. 

The result of the calculation of R s = V s /I s is shown in Fig.|9| for the following 
parameters : I s = 1/iA, H c2ab {0) = 674T |TJ, £ ab (0) = lAnm Q, k = M£ = 58 
fTTf, n = 4.5 x 10 6 , T c = 8AAK, L = 0Amm,A = (2nm) x (lnm), 7 = 0.067 and 
n s = 8. According to Fig.|l|, the peak does not contribute to p(T) at temperatures 
above approximately 84 K, which corresponds to V s = 0. The LA-MH theory does 
not apply at temperatures very close to T c because of the condition for the applied 
current I s which must be smaller than I c (Eq.(^|)). Therefore in the calculation we 
used T c about 0.4 K higher. Good agreement between the experimental data and the 
calculated resistivity versus temperature was obtained (see Fig.|S]). 

The experimental data show that the reduction of the resistive peak magnitude 
in p(T) and the width of the peak in p(6) at 9 = 0° occurs when the magnetic 
field direction is rotated from the ab-planes (8 = 0°) towards the c-axis (8 = 90°) (see 
Fig.|I|). When B is rotated from 8 = 0° position, its c-axis component B sin 8 increases 
and the critical current I c in the ab-planes decreases. The resistive peaks in p(T) and 
p(8) at 8 = 0° also decrease in magnitude with an increasing applied transport current 
I s (see Fig.^). We verified, using numerical calculations that according to the LA-MH 
theory, the phase-slip voltage decreases with an increasing I s and decreasing I c in the 
limit of very small currents (see Fig.|T0|). 
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5 SUMMARY AND CONCLUSIONS 



We investigated the resistive peak anomaly in underdoped YBCO, which was observed 
in both the temperature dependence of resistivity p(T) and the angular dependence 
of resistivity p(9) in an applied magnetic field B. The resistive peak anomaly in p(T) 
decreases with an increasing B (applied parallel to the c-axis) and with an increasing 
applied transport current I s . On the other hand, the width of the resistive peak in 
p(8) at 9 = 0° decreases with an increasing B, and its magnitude decreases with 
an increasing I s . The resistive peak anomaly in p(T) and its dependence on B and 
I s show striking qualitative similarities to those exhibited by LTSC wires, and some 
LTSC thin films and HTSC crystals. The YBCO film that we analyzed has resistivity 
much lower than YBCO crystals studied by Mosqueira et al [§], suggesting that the 
resistive peak anomaly is not directly related to the absolute value of the normal state 
resistivity. The anomaly can not be explained by the c-axis misalignments, since they 
would eliminate the sharp minimum in p(9) at 9 = 0° (B || ab-planes) observed in 
regions I and III. 

We analyzed the data in terms of three different models that were developed 
in the past to explain the resistive anomalies. The 2D-resistor model and the flux 
motion models are inadequate to explain fully our data including the dependence on 
a magnetic field and an applied transport current. The phase-slip (LA-MH) model 
provides the best qualitative and quantitative description of the observed resistive 
anomalies and their behavior as a function of temperature T, magnetic field B, the 
angle between B and the ab-planes, and the applied transport current I s . This model 
can be applied under the assumption that the current flows through ID filaments. The 
assumption about the filamentary flow of the current was supported by the data of 



Browning et al fL2 |, which revealed that in spite of a very large variation of the oxygen 
content in YBCO single crystals, they have very sharp superconducting transitions, 
very high T c , and low resistivity. 
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Figure 9: (a) The assumed normal state resistivity without phase-slip resistivity 
contribution in a zero magnetic field; (b) The calculated phase-slip resistivity ;(c) 
The fitting to the experimental data obtained by superposition of the normal state 
resistivity in (a) and the calculated phase-slip resistivity in (b). 
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Figure 10: The calculated phase-slip resistivity as a function of the normalized 
critical current (I c /h) at different applied currents I s . 



27 



